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Abstract: A common method for measuring the electromagnetic properties of super- 
conductors is to measure their complex magnetic ac susceptibility x as a function of 
frequency u and amplitude H , and of temperature T and applied dc magnetic field H^ c 
as parameters. The basic theory of the linear x{ u ) an d nonlinear x(Hq,u) is outlined 
for various geometries, e.g. disks, rings, and strips of thin films or thicker platelets in 
a perpendicular magnetic field. It is shown how xO^O explicitly depends on the linear 
resistivity p ac (uj) = E/J or penetration depth X ac (u), and x(H ,uj) on the nonlinear 
current- volt age law E(J, B) where E(r), J(r), B(r) are the local electric field, current 
density, and induction. The dependence of E(J, B) on T and on various material proper- 
ties like pinning forces or pinning energies, structural defects and granularity, leads to an 
implicit dependence of x on these parameters. 

£ ■ 1. Introduction 

q \ The magnetic response of type-II superconductors may be considered under various 



>< 



aspects. First, one may be interested in the microscopic motion of the Abrikosov vor- 
tices, i.e. in their pinning, thermally activated depinning, crossing of surface barriers, and 



mutual annihilation. Second, within continuum theory these microscopic processes are 
described by a current-voltage characteristic E = E(J, B) (the local electric field caused 
by the current density J and depending also on the local induction B) and by a reversible 
magnetization curve H = H(B), which in sufficiently large B may be approximated by 
H = B//xq, but at smaller B may lead to a geometric barrier for flux penetration (Zel- 
dov et al. 1994, Brandt 1999). In general, these material laws are anisotropic when the 
material is anisotropic, when anisotropic defects are introduced, or when the Hall effect 
is taken into account. In the isotropic case one has H = H(B)~H/H and (in simple ge- 
ometries) E = p(J, B)3 with p = E/J. Third, the geometry of the experiment, i.e. the 
specimen shape and orientation of the applied magnetic field and/or transport current 
have to be accounted for. The geometry determines the local profiles B(r) and J(r) and 
the global response, i.e. the magnetic moment and voltage drop. In some cases the ge- 
ometry modifies the appropriate law E(J,B). For example, in isotropic films the critical 
current density J c becomes anisotropic when a strong magnetic field is applied parallel to 
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the film (Schuster et al. 1997). Fourth, one may consider the vortex state as a "phase" 
which undergoes transitions from a strong pinning state ("vortex glass") to a thermally 
depinned state when in the BT plane an "irreversibility line" T rev (B) or "glass line" T g (B) 
is crossed (Blatter et al. 1993, Brandt 1995), or which changes from solid to "liquid" at 
a "melting line" T m (B) where the flux-line lattice melts (Zeldov et al. 1995, Majer et 
al. 1995, Schilling et al. 1996, Welp et al. 1996, Roulin et al. 1998, Sasagawa et al. 
1998, Dodgson et al. 1998). In such a "universal" description the detailed microscopic 
behavior of the vortices remains less clear; in particular, the microscopic interpretation of 
the length and time scales which diverge at the glass temperature T g (B) is still open. 

At low temperatures T, the current-voltage curve E(J, B) is highly nonlinear, typically 
E = E c ■ (J/J c ) n with large exponent n(B,T) ^> 1, and one arrives at the usual Bean 
model (Bean 1964, Campbell and Evetts 1972) with critical current density J C (B,T). 
With increasing T, the exponent n decreases and one observes flux creep, i.e. nonlinear 
relaxation of currents and magnetic moment, see Sec. 2. At high T above the irreversibil- 
ity line, thermal depinning leads to a linear E(J) which in general is frequency dependent 
and complex. An interesting task is the detailed measurement of E(J, B) of thin super- 
conductor films at various temperatures. A nontrivial reversible H(B) with a jump at 
the lower critical field H c \ (calculated, e.g. from the Ginzburg-Landau equations) leads to 
surface barriers, which may be microsopic as predicted by Bean and Livingston (1964), 
or geometric, caused by the non-ellipsoidal shape of the superconductor. 

Some geometry problems have been solved exactly, e.g. the Bean model for long su- 
perconductors in parallel field (Bean 1964) and for thin disks (Mikheenko and Kuzovlev 
1993) and long strips (Brandt et al. 1993) in perpendicular field, and in strips with trans- 
port current without (Norris 1970) and with (Brandt and Indenbom 1993, McElfresh et 
al. 1994) applied magnetic field. Known are further the linear (Clem et al. 1976, Kes 
et al. 1989) and nonlinear (Rhyner 1993, Gilchrist and Dombre 1994) relaxation (or flux 
diffusion) in longitudinal geometry where demagnetization effects are absent, and the lin- 
ear and nonlinear reversible response of homogeneous ellipsoids, which is obtained from 
the longitudinal results by introducing a demagnetization factor. The magnetic response 
of thin strips and disks (Brandt 1994a,b) in a perpendicular field is easily computed by 
solving a one-dimensional (ID) integral equation for the current density, while thin rect- 
angular (or arbitrarily shaped) films and thick strips and disks are two-dimensional (2D) 
problems, which still can be computed on a PC (Brandt 1995a, 1996b, 1998). 

Recently two different algorithms were presented (Labusch and Doyle 1997, Doyle et 
al. 1997, Brandt 1999) which account for both constitutive laws E = E(J, B, r) and 
H = H(B, r) and which apply to all geometries, i.e. to superconductors of any shape. 
These universal methods allow to compute the geometric "edge barrier", (Zeldov et al. 
1994), e.g. in strips and circular disks (or cylinders) of finite thickness in perpendicular 
field. For the reversible magnetization curve H(B) any model may be used which ap- 
proximates the exact Ginzburg-Landau result (Brandt 1997a). If one is interested in the 
behavior not too close to the upper critical field H c2 , one may use a simple hyperbola: 



H(B) rs [H%i + (-B//xo) 2 ] 1//2 sign(_B), for B > jjlqH c i and with undefined reversible field 
-H cl < H < H cl at B = 0. 

As mentioned above, really 3D problems (e.g. the superconducting cube) have not been 
computed so far. The difficulty is not only of numerical nature, but it is even not clear 
which material laws E(J, B) are reasonable when J is not parallel to B. One may introduce 
two critical current densities J c ± and J c \\ for currents _L B and \\B. In general, at low 
B <C B C 2 (the upper critical field) one should have E oc B since each vortex contributes 
independently to the resistivity. 

Strictly spoken, even the flux penetration into thin films is a 3D problem if the specimen 
is not of exactly circular shape or is not an infinitely long strip (Mikitik and Brandt 
1999). This is so since during the penetration of flux, the orientation of the sheet current 
3 s (x,y) = J 3(x,y,z) dz changes with time. Therefore, the orientation of the vortices 
which penetrate from the flat surfaces also changes, since these vortices are perpendicular 
to 3 S . This means that at any point x, y away from the symmetry axes and from the 
specimen edge or center, in principle the vortex orientation may be a function of the 
depth z, i.e. the vortex arrangement may be twisted. This rotation is most pronounced 
at moderate magnetic fields; it is absent when the applied field is above the field of full 
penetration, or when the film thickness d is smaller than the magnetic penetration depth 
A, since in these cases the vortices are nearly along z. 

2. Definition of AC Susceptibilities 

The global magnetic response of a superconductor of arbitrary shape in a time dependent 
homogeneous magnetic field H a (t) may be characterized by its (dipolar) magnetic moment 



m(£) = - / r x J(r, t) d 3 r , 



were J(r, t) is the current density inside the specimen of volume V. Extensions to inhomo- 
geneous H a (r, t) and to higher (e.g. quadrupolar) moments usually are not required when 
an experiment is properly designed to obtain the material laws of the specimen. These 
material laws may be linear or nonlinear. In the linear case, the conductor is completely 
described by its linear resistivity p ac {uj) = iuj^\ 2 ac) which in general is complex and fre- 
quency dependent, or equivalently, by a conductivity o" ac (cj) = 1/ p ac {uj) or a complex 
penetration depth A ac (cj) = (pac/itopo) 1 ^ 2 . For example, an applied field H a (i) = H e lu}t 
generates a local current density J(r, t) = 3o{r)e tuJt and an electric field E(r, t) = Eo(r)e lu}t 
where Jo(r) and Eo(r) = p a cJo( r ) are complex amplitudes. The magnetic field then 
penetrates into a conducting half space x > as H(x,t) = H a (t) exp[—x/ \ ac (u)]. For 
magnetic materials with linear permeability fi(u) = B/(fi H), in these formulae /xo has 
to be replaced by fi(u)fiQ. A superconductor may exhibit such a linear response above a 
temperature where thermal depinning of vortices, or thermally assisted flux flow (TAFF) 
occurs (Kes et al. 1989). 

If pinning is strong, or if the reversible magnetization curve H = H(B) with a finite lower 
critical field H c i is accounted for, the magnetic response is nonlinear. At sufficiently large 



induction B ^> HqH c \ one may put B = /ioH. In this case pinning and flux creep may be 
modeled by a nonlinear resistivity law E = p(J)J, e.g. by the model p(J) = Po(^/^c) n_1 
with critical current density J C (B) and creep exponent n(B) > 1, as mentioned above. 
This model results when in the electric field E(J) = E exp[—U(J)/kT] an activation 
energy of the form U{ J) = Uq ln( J c / J) is inserted, yielding E = E c (J/J c ) n with n = 
U /kT. 

Interestingly, when E oc J n then during creep in the fully penetrated state, after some 
transient time ti the electric field exhibits universal behavior, i.e. the r and t dependences 
separate: E(r,t) oc /(r)/(ti + t)"/( n_1 ) m /(r)/i, where the profile /(r) depends only 
on the specimen shape (Gurevich and Brandt 1994, Brandt 1996a). As a consequence, 
J(r, t) oc /(r) 1 / n /(t 1 + t) 1 /^-- 1 ) oc const — lnt if n ^> 1, and thus the magnetic moment 
relaxes as m oc const — In t. During flux creep the magnetic response to a small ac magnetic 
field is linear, though the underlying constitutive law E oc J n is highly nonlinear (Gurevich 
and Brandt 1997, Brandt and Gurevich 1996). 

From the calculated current density J(r,H a ) the magnetic moment m(H a ) (1) is ob- 
tained. If H a (t) is cycled sinusoidally, H a (t) = H sin(ut), then m(t) (the component of 
m along H a ) performs a hysteresis loop from which nonlinear complex ac susceptibilities 
X^(H , u>) = x'/j, ~ ixl (/■* = 1, 2, 3, . . .) may be defined as 

X ,(H , u) = -^=- r m (t) e-^ d{oot) . (2) 

TT-tlo JO 

For small amplitudes H — > this yields xi(Hq,w) = to'(O) = \im Ha _> dm(H a )/dH a . A 
convenient normalization of the Xn(Ho,co) is thus Xn ~ > X /j /|'?tz. / (0) | , which at H — ► 
yields the ideal diamagnetic susceptibility Xi(0,u;) = —1 if m'(0) < 0. In the critical state 
model, with general J C (B) or with Bean's J c =const, the susceptibilities x M (if , w) do not 
depend on cu, i.e. all cycling rates give the same loop m(H a ) and same Xn(Ho). 

In the opposite limit of linear response, no higher harmonics are generated, i.e. Xn = 
for /x ^ 1, and %i = x(^) does not depend on the amplitude. In complex notation with 
H a (t) = Hoe luJt (only the real parts of H a , H, B, E, J, and m have physical meaning) the 
linear ac susceptibility x(w) — x' ~ ^x" is defined as 

x(w) = 4- r^) ^ d ^) ■ ^ 

In the general nonlinear case the x^{Hq^oo) depend on both Hq and uo. In the use- 
ful model E(J) = E c (J/J c ) n , the x^Hq^oj) depend only on combinations of the form 
Hq/uj 1 /^ 1-1 ^ or lu/Hq^ 1 or any function of these ratios. Thus, x(Ho,u) for different fre- 
quencies uj is obtained by rescaling the amplitude axis, x(Ho,cu) = x{Ho/c 1 ^ n ~ 1 \oj) 
for any constant c. This scaling to a good approximation applies also to other E(J) 
laws if these are sufficiently nonlinear and if the effective exponent n is defined as n = 
<9(ln E)/d(ln J) taken at J = J c where J c is the typical current density of the experiment. 
The nonlinear susceptibility thus depends only on one variable combining amplitude and 
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frequency, and further on an effective exponent n and on the geometry. 

3. Some Linear AC Susceptibilities 

The linear susceptibilities x{°°) = x'~^x" or permeabilities p{uo) = x(cj) + 1 of conductors 
with arbitrary shape and arbitrary complex ac resistivity p ac or penetration depth X ac = 
\p ac {uj)/iujjiQ\ 1 / 2 in a homogeneous magnetic field H a {t) = H§e wt may be obtained by 
solving a linear diffusion equation for J(r,t) or H(r,t) with diffusivity D = p ac j 'p . In 
this way one finds for infinite slabs (width 2a) in a parallel field (Kes et al. 1989) and 
long cylinders (Clem et al. 1976) or spheres of radius a (London 1961) in an axial field: 

(4) 

XcylM = SS " '' (5) 

, 3 coth u 3 

XsphercM = 2 ~ 1 ' W 

where w = a/\ ac = [iujpoa 2 / ' p{uo)} 1 / 2 is complex and the definition x( u ) = Tn(u)/\m{u — ► 
oo) | was used. Jo and Ji are modified Bessel functions. At high frequencies (\u\ ^> 1) one 
has for /i = x + 1: A^siab ~ l/«, At cy i ~ 2/u, and yU sphcrc w 3/u. 

Interestingly, long cylinders in perpendicular H a (t) yield the same normalized suscepti- 
bility Xcyi(^) (5) as cylinders in parallel field, but the magnetic moment is twice as large: 
m cyiW = 2^c y i = 27ra 2 i^ a (t)x cy i(cij). Since this response is linear, one may superimpose 
both solutions and thus finds the magnetic moment for long cylinders in an applied field 
H a (t) inclined at an arbitrary angle 9 away from the cylinder axis, m cy i = (mjL, w^L), 

m cy i(£;0) = 7ra 2 iJ a (t)xcyi(^)(cos6', 2sin6>). (7) 

The current density in the cylinder is J = (J v , J\\), 

J(r,<p,t;0) = ^^{^h cose, 2sinflcosy9). (8) 

In general, the linear susceptibility for any geometry may be written as an infinite sum, 
or approximated by a finite sum, of the form 



V 



Here A u {v — 1, 2, . . .) are the eigenvalues of an eigenvalue problem, the b v ("dipole 
moments") are integrals over the eigenfunctions /^(r), and the complex variable w is pro- 
portional to iuj J ' p ac (uS). The sum in the denominator of Eq. (9) provides the normalization 
x(oo —> oo) = —1 (ideal diamagnetic screening). 



For thin film disks and strips with diameter or width la and thickness 2b <C la in 
perpendicular field the A u and b u are tabulated by Kotzler et al. 1994 and Brandt 1994c, 
and the main variable is 

iup ab ab 
w = — = — — - = iwt(w) . 10 

-np ac {uj) it\i c 

Tables of the A u and b v for thick disks or cylinders of arbitrary length in axial field are 
given by Brandt 1998. For infinite bars with rectangular cross-section \x\ < a, \y\ < b, in 
perpendicular H a (t) along y, the variable w (10) applies also and the eigenvalue problem 
reads (Brandt 1996b) 

A 



Uv) = -^fdx' Jdy'\n\v-v'\f v {v'), (11) 



where r = (x,y), r' = (x',y f ). With the normalization / f^f u d 2 r = 5^ the "oscillator 
strengths" are in this case 



h .xfv(x,y)dr. (12) 

For practical purposes, a finite number of terms v — 1 . . . N in the sum (9) is sufficient. 
When the (real and positive) numbers A u and b v are known for a given geometry, then x{ u ) 
(9) may be calcuated for any complex resistivity p ac (u>). By inverting this relationship 
between the two complex functions x(w) and p ac (uj) numerically, the complex resistivity 
p ac (oo) may be obtained from measured ac susceptibilities as done by Kotzler et al. 1994. 

The linear x'{ u ) an d x"^) °f diks or cylinders with Ohmic resistivity p a c(^) — P =const 
in axial H a (t) are shown in Fig. 1 for various aspect ratios b/a = ... 10 (a = radius, 
lb = length) versus the reduced frequency uor with r = p ab/(7rp). Note that in this 
double logarithmic plot for sufficiently thin disks {b <^ a) the p(uj) = xi 00 ) + 1 with 
increasing frequency crosses over from the behavior in perpendicular geometry (nonlocal 
flux diffusion) (Brandt 1994c), 

p! — x' + 1 oa 1/cj, p" = x" oc ln(const • ur)/u , (13) 

to the parallel geometry (local flux diffusion), 

H'pafjf'ocl/y/u. (14) 

At all aspect ratios b/a > the real and imaginary parts at large frequencies coincide, 
p'(uj) = p"(oo). The physical reason for this finding is that above the frequency where 
the skin depth 5 = V2\ ac — (Zp/i^Po) 1 ^ 2 coincides with the half thickness 6, the Ohmic 
conductor nearly behaves like an ideal diamagnet (or superconductor in the Meissner 
state), screening almost all magnetic flux from the interior of the conductor. The magnetic 
field lines thus have to flow around the bar or cylinder such that B is parallel to the 
specimen surface everywhere. Therefore, at high frequencies, thin (and thick) Ohmic 
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(and non-Ohmic) conductors in a perpendicular magnetic ac field behave as if the field 
were applied parallel. 

4. Nonlinear Magnetic Response of Thin Films 

The nonlinear susceptibility of thin superconductor films with given material law in a 
perpendicular magnetic field, in general may be obtained only numerically, e.g. by time- 
integrating an integral equation for the current density J(r,t), inserting this into Eq. (1) 
for the magnetic moment m(t), and then evaluating the Fourier integral (2) for x{Hq,lo). 
Such results are presented for thin film disks and rings by Brandt 1997b and for thick 
disks (or short cylinders) by Brandt 1998. 

For narrow rings or thin-walled hollow cylinders, analytical expressions for x{Hq) are 
available in the Bean limit (Brandt 1997b). In this case (and also when a ring contains 
one or several weak links) the virgin magnetization curve is a straight line, m(H a ) = 
—H a m sa t/H p (ideal screening or Meissner state) which saturates to a horizontal line, 
l m l — "^sat at \H a \ > H p , the field of first penetration of flux into the hole. This 
is so since the current in the ring is limited to a critical value I c . For a narrow ring 
with width w = a — a\ <C R = (a + a±)/2 (a, a,\ are the outer and inner radius) and 
thickness d <C R one has I c = J c wd, m sat = ttR 2 I c , H p rs I c ln(5R/x)/irR where the 
inductivity L ps fi Rln(5R/w) was used (Brandt 1997b). For such rings or tubes, in cycled 
H a (t) = Hosinuit the magnetization loop at small amplitudes Hq < H p degenerates to a 
straight line, yielding a constant and real x{Hq) = —m sa , t /H p , while at large Hq > H p 
the loop m(H a ) is a parallelogramme with height 2m sat . This yields the normalized ac 
susceptibility x(h) = x(H)H p /m sat = x' — ix", with h = H /H p and s = 2/h — 1, 

X\h) = -1, 

x "(h) = 0, h<l, 

X'(h)= ~\-\ arcsin s-\syj\- s 2 , 

Ah -I 1 - s 2 

X'\h)= lf V = ^A h>l, (15) 

see the curve a\ja = 1 in Fig. 2. The polar plot {x" versus x' with h as parameter) of 
the ring susceptibility (15) is symmetric, i.e. x"{x') yields the same curve as x"{~ 1 ~x')- 
The maximum of Xmax = V 71 " = 0.318 occurs at h — 2 (at s — 0). For large amplitudes 
h = H /H p > 1 one has x(h) w -1.69/h 3 ^ 2 and X "(h) w 4/(irh). 

Very good approximate expressions are available within the Bean critical state model 
for thin strips (width 2a, thickness d <^C a, length L ^> a), disks (radius a), and ellipses 
(semi-axes a and b, excentricity e = b/a < 1) (Mikitik and Brandt 1999). In these three 
cases and also for square and rectangular films, the virgin magnetization curve with less 
than 1% error is 

m(H a ) = -TO sat tanh(if a /ifi) , (16) 



with H\ = TOgat/ |m'(0) , and explicitly, 










m sat = J c dLa 2 , 




Hi = J c d/7T 
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^sat — —JcdO, , 
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The formulae (19) for ellipses with e = 


= b/a 


< 1 reproduce the results for disks (e 


-1, 



radius a = b) and strip (e — > 0, length 2a, width J8/3b). From the virgin curve m(H a ) 
(16), the Bean magnetization loops with amplitude Hq follow as 

mi (H a , H ) = m(H ) + 2m( Ha ~ ^ ) (20) 



(decreasing H a ), and m-[(H a , H ) = —mi(—H a ,H ) (increasing H a ). Inserting this into 
definition (2) one obtains the universal nonlinear ac susceptibility x(h) of films normalized 
to x(0) = —1 and depending on h = Hq/Hi: 



oj /"t/2 r , ~\ — sin<j3\ i 

X(h) = —r [2 tanh (h -) - tanh h] e~ lip dip 



2i f n / 2 r , /, 1 — sin<z>> 
2tanh(/? 

-jr/2 

From this universal function the not normalized x{H$) is obtained as 



(21) 



X(tf ) = ^x(f). (22) 

Figure 2 shows some normalized ac susceptibilities x(h) for thin film rings and for the 
disk and strip in the Bean limit, i.e. for large creep exponent n — > oo. Note that the 
curves for disks and strips practically coincide. 

The nonlinear ac susceptibility of disks with various creep exponents n > 3 is depicted 
in Fig. 3. Here E(J) = E c (J/J c ) n was assumed, the circular frequency was uo = ooi = 
2E c /(n J c da), and the amplitude H is in units J c d (d = tickness, a = radius of the 
disk). The plots in Fig. 3 apply to any frequency u if the unit of Hq is changed to 

(u/uifl^-^JJ. 

5. Final Remarks 

The measured ac susceptibilities depend on the temperature T indirectly via the T 
dependence of the material properties. The nonlinear x{Hq,uj) depends on J C (B,T) and 
on the creep exponent n(B,T) = Uo/kT, see above. The linear x{°°) depends on the 
variable w, Eq. (10), which contains the linear complex ac resistivity p ac {^, B, T), see e.g. 
Coffey and Clem 1990 and Brandt 1990 for theories of p ac (^) of type-II superconductors. 
The linear xi 00 ) usually is measured in a bias dc magnetic field H^ c ^> H , but the 
nonlinear x{Hq-,oo) discussed above assumed the absence of a bias field and has to be 
recomputed if H^ c > 0. 



Furthermore, if H a (t) is too small, the effects of a finite H c \ should be considered, which 
were disregarded here but may lead to a geometric barrier, in particular in thin films, 
as discussed in the introduction. The edge barrier may be suppressed by using small 
measuring coils close to the film, see e.g. Gilchrist and Brandt 1996 for a detailed theory 
and analytic formulae. The numerical program described by Brandt 1999 in principle 
accounts for arbitrary reversible magnetization H(B) and thus for H rA and for the edge 
barrier in strips or disks of constant thickness. It may also be generalized to allow for dc 
and ac transport currents in strips. 

The problem of ac losses in thin strip superconductors with transport current is related 
to the problem of ac susceptibilities. It is conceptually difficult since with increasing ac 
amplitude the losses may be caused first by flux jumping over the edge barrier, then by 
bulk pinning (Bean critical state), then (after full penetration) by the rapidly increasing 
electric field E oc J n . The hysteresis losses caused by the edge barrier in superconductor 
strips without bulk pinning were recently calcuated by Clem and Benkraouda 1998. 

5. References 



Bean, C. P., 1964, Rev. Mod. Phys., 36, 31. 

Bean, C. P., and Livingston, J. D., 1964, Phys. Rev. Lett, 12, 14. 

Blatter G., Feigel'man M. V, Geshkenbein, V. B., Larkin, A. I., and Vinokur, V. M. 
1994, Rev. Mod. Phys., 66 1125. 

Brandt, E. H., 1991, Phys. Rev. Lett, 66, 3213. 

., 1994a, Phys. Rev. B, 49, 9024. 

., 1994b, Phys. Rev. B, 50, 4034. 

., 1994c, Phys. Rev. B, 50, 13833. 

., 1995a, Phys. Rev. B, 52, 15442. 

., 1995b, Rep. Prog. Phys., 58, 1465. 

., 1996a, Phys. Rev. Lett, 76, 4030. 

., 1996b, Phys. Rev. B, 54, 4246. 

., 1997a, Phys. Rev. Lett, 78, 2208. 

., 1997b, Phys. Rev. B, 55, 14513. 

., 1998, Phys. Rev. B, 58, 6506, 6523. 

., 1999, Phys. Rev. B, 59, 3369 

., and Gurevich, A., 1997, Phys. Rev. Lett, 76, 1723. 

., and Indenbom, M., 1993, Phys. Rev. B, 48, 12893. 

., Indenbom, M., and Forkl, A., 1993, Europhys. Lett., 22, 735. 
Campbell, A. M., and Evetts, J. E., 1972, Adv. Phys., 72, 199. 



Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 


Brandt 


E. 


H.. 



Clem, J. R., and Benkraouda, M., 1998, Phys. Rev. B, 58, 15103. 

Clem, J. R., Kerchner, H. R., and Sekula, T. S., 1976, Phys. Rev. B, 14, 1893. 

Coffey, M., and Clem, J. R., 1991, Phys. Rev. Lett, 67, 386. 

M. J. W. Dodgson, M. J. W., Geshkenbein, V. B., Nordborg, H., and G. Blatter, 1998, 

Phys. Rev. Lett, 80, 837; Phys. Rev. B, 57, 14498. 
Doyle, T. B., Labusch, R., and Doyle, R. A., 1997, Physica C, 290, 148. 
Gilchrist, J., and Brandt, E. H., 1996, Phys. Rev. B, 54, 3530. 
Gilchrist, J., and Dombre, T., 1994, Phys. Rev. B, 49, 1466. 
Gurevich, A., and Brandt, E. H., 1994, Phys. Rev. Lett, 73, 178. 
Gurevich, A., and Brandt, E. H., 1997, Phys. Rev. B, 55, 12706. 

Kes, P. H., Aarts, J., van den Berg, J., van der Beek, C. J., Mydosh, J. A., 1989, 
Supercond. Sci. Technol, 1, 242. 

Kotzler, J., Nakielski, C, Baumann, M., Behr, R., Goerke, F., and Brandt, E. H., 1994, 

Phys. Rev. B, 50, 3384. 
Labusch, R., and Doyle, T. B., 1997, Physica C, 290, 143. 

London, F., 1961, Superfluids II. Macroscopic Theory of Superconductivity (New York: 
Dover), p. 35. 

Majer, D., Zeldov, E., and Konczykowski, M., 1995, Phys. Rev. Lett., 75, 1166. 

McElfresh, M., Zeldov, E., Clem, J. R., Darwin, M., Deak, J., and Hou, L., 1995, Phys. 

Rev. B, 51, 9111. 
Mikheenko, R N., and Kuzovlev, Yu. E., 1993, Physica C, 204, 229. 
Mikitik, G. P., and Brandt, E. H., 1999, Phys. Rev. B (submitted). 
Norris, W. T., 1970, J. Phys. D: Appl. Phys., 3, 489. 
Rhyner, J., 1993, Physica C, 212, 292. 
Roulin, M., Junod, A., and Walker, E., 1998, Phys. Rev. Lett, 80, 1722. 

Sasagawa, T., Kishio, K., Togawa, Y., Shimoyama, J., and Kitazawa, K., 1998, Phys. 

Rev. Lett, 80, 4297. 
Schilling, A., Fisher, R. A., Phillips, N. E., Welp, U., Dasgupta, D., Kwok, W. K, and 

Crabtree, G. W., 1996, Nature, 382, 791. 
Schuster, Th., Kuhn, H., Brandt, E. H., and Klaumiinzer, S., 1997, Phys. Rev. B 56, 

3413. 
Welp, U., Fendrich, J. A., Kwok, W. K., Crabtree, G. W., and Veal, B. W., 1996, Phys. 

Rev. Lett, 76, 4809. 
Zeldov, E., Larkin, A. I., Geshkenbein, V. B., Konczykowski, M., Majer, D., Khaykovich, 

B., Vinokur, V. M., and Shtrikman, H., 1994, Phys. Rev. Lett, 73, 1428. 

Zeldov, E., Majer, D., Konczykowski, M., Geshkenbein, V. B., Vinokur, V. M., and 
Shtrikman, H., 1995, Nature, 375, 373. 



10 



+ 0.1 

"x 

II 

0.01 



X 

II 



0.1 



0.01 




/ / 

J_LI_l£ I I I 



0.1 



10 
COT 



100 



1000 



Figure 1: The real (top) and imaginary (bottom) parts of the linear ac permeability 
yU = X+l=/i' — ifi" of Ohmic cylinders with rectangular cross-section 2a x lb in an axial 
magnetic ac field for aspect ratios b/a = (thin disk), b/a = 0.01, 0.03, 0.1, 0.3, 1, 2, 
and 3 (short cylinders) computed from the sum (9) (solid lines). The dashed curves give 
the analytic expression (5) for long cylinders with b/a = 3, 5, and 10, with the variable 
u = a/X ac = {iuoTTxa/b) 1 / 2 inserted, where r = fioab/(irp) is the time scale. For b/a = 3 
the numerical and analytical curves are very close. To facilitate comparison, some yJ 
curves are repeated as dotted lines in the lower plot. 
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Figure 2: Top: Real and imaginary parts of the nonlinear complex susceptibility x{Hq) = 
x' ~ i>x" °f thin strips, disks, and rings with various ratios a^ja of the inner and outer 
radius, plotted versus the amplitude H of the ac field. Bean model with constant J c . The 
unit of H is Hi = J c d/n for the strip and Hi = 7iJ c d/8 for the disk. For rings H is in 
units of a penetration field H' pa H p which gives best fit to the sharp rise in x" °f the ideal 
ring: H' p /J c d = 0.126 (0.210) for oi/o = 0.9 (0.8), while from Table 1 in Brandt 1997b 
one has H p /J c d = 0.136 (0.239). For the narrow ring with width w = a - ai -C a one has 
H' = H p , H p /J c d pa (w/tto) ln(5a/w). Bottom: Polar plots of the same X = x' ~ ^x" ■ For 
the narrow ring x' an d x" [Eq. (15)] are discontinuous at Hq = H p and the polar plot is 
symmetric. 
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Figure 3: Top: Real and imaginary parts of the nonlinear complex susceptibility x{Hq) — 
x' ~ ix" °f a thin disk plotted versus the ac amplitude Ho for various creep exponents 
n = 3, 5, 11, 51, and 201 in the law E = E c (J/J c ) n with constant J c . Here Hq is in units 
J c d. Bottom: Polar plots of the same x = x' ~~ i>x" ■ The dashed line shows the Bean 
model (n — > cxd) for the disk. 
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